Excited-state phase transition leading to symmetry-breaking steady states in the 

Dicke model 
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We study the phase diagram of the Dicke model in terms of the excitation energy and the radiation- 
matter coupling constant A. Below a certain critical value A c of the coupling constant A all the energy 
levels have a well-defined parity. For A > A c the energy spectrum exhibits two different phases 
separated by a critical energy E c that proves to be independent of A. In the upper phase, the energy 
levels have also a well defined parity, but below E c the energy levels are doubly degenerated. We show 
that the long-time behavior of appropriate parity-breaking observables distinguishes between the 
different phases of the energy spectrum of the Dicke model. Steady states reached from symmetry- 
breaking initial conditions restore the symmetry only if their expected energies are above the critical. 
This fact makes it possible to experimentally explore the complete phase diagram of the excitation 
spectrum of the Dicke model. 
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Introduction.- The rapid development of experimen- 
tal techniques controlling ultra-cold atoms has given 
rise to a great breakthrough in the physics of quan- 
tum many-body systems. A logical outcome has been 
the increase of the interest in certain phenomena, such 
as non-equilibrium dynamics and quantum phase transi- 
tions (QPTs). Also, it has entailed the revival of well- 
known physical models as that formulated by R. Dicke, 
which describes the interaction between an ensemble of 
two-level atoms with a single electromagnetic field mode, 
as a function of the radiation- matter coupling Its 
most representative features are a second order QPT, 
which leads the system from a normal to a superradiant 
phase, characterized by a macroscopic population of the 
upper atomic level Q , the emergence of quantum chaos 
and the spontaneous symmetry breaking [j, Q ■ Although 
this model has been extensively studied from many points 
of view, there still exists a heated controversy about its 
significance in real physical systems. A no-go theorem 
was formulated in the seventies, stating that the super- 
radiant phase transition cannot occur in a general sys- 
tem of atoms or molecules interacting with a finite num- 
ber of radiation modes in the dipole approximation 
In addition, it is not clear if this theorem also forbids 
the superradiant transition in other realizations of the 
Dicke model, like in circuit QED [6j. On the contrary, 
this transition has been experimentally observed with a 
superfluid gas in an optical cavity, giving rise to a self- 
organized phase [7] . All these facts have turned the Dicke 
model into a multidisciplinary hot topic, involving differ- 
ent branches of physics. As a consequence, there exists 
an intense theoretical research; a few representative ex- 
amples concern non-equilibrium QPTs [8j , thermal phase 
transitions in the ultrastrong-coupling limit or equi- 
libration and macroscopic quantum fluctuations fioj ]. 

In this Letter we explore the phase diagram of the 



Dicke model as a function of two control parameters: the 
radiation-matter coupling constant A and the energy E 
of its eigenstates. We show that the energy spectrum can 
be divided into three different sectors or phases separated 
by certain critical values A c and E c . For A < A c we find 
that parity is a well defined quantum number at any exci- 
tation energy. The situation is rather different if A > A c . 
Below a certain critical energy E c all the energy levels 
of the system are doubly degenerated, and, as a conse- 
quence, the parity symmetry of each level can be broken. 
Above the critical energy E c there are no such degen- 
eracies and parity is again a good quantum number. To 
some extent we can say that beyond A c the excited energy 
levels up to energy E c inherit the properties of the super- 
radiant phase, characteristic of the ground state. We also 
show that this phase diagram entails measurable effects 
in the long-time dynamics of certain observables. In- 
deed, if one prepares the system in a symmetry-breaking 
ground state of the superradiant phase and performs a 
quench to a non-equilibrium state, then the symmetry 
of the final steady state remains broken if its energy is 
below E c , whilst it is restored in the opposite case. As a 
consequence, parity non-conserving observables relax to 
steady values different from zero only if the energy of the 
non-equilibrium state is below the critical. This fact con- 
stitutes an unheralded characteristic of the Dicke model, 
that can be accessible to experiments and shed some light 
over the current controversy about the relevance of the 
critical behavior of this model in real physical systems. 

The Dicke model. - The Dicke Hamiltonian can be writ- 
ten as follows 

2A 

: J z + uo^a H — j= (a i + a) J x , (1) 
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where and a are the usual creation and annihilation 
operators of photons, J = ( J x , J y , J z ) is the angular mo- 
mentum, with a pseudo-spin length J = N/2, and N 
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is the number of atoms. The frequency of the cavity 
mode is represented by u> and the transition frequency 
by uj . Finally, the parameter A is the radiation-matter 
coupling. Throughout all this Letter, we take H=l, and 
w = w = l. The parity II = e H J + J *+^ a ) i s a conserved 
quantity, due to the invariance of H under J x — » —J x and 
a — » —a 0, and thus all the eigenstates are labeled with 
positive or negative parities. The system undergoes a 
second-order QPT at A c = which separates the 

so called normal phase (A < A c ) from the superradiant 
phase (A > A c ) [2j . In the latter the ground state becomes 
doubly degenerated and parity can be spontaneously bro- 
ken — because of the fluctuations, the system can evolve 
into one particular ground state without a well-defined 
parity [4j. 

It is important to note that the system has a finite 
number of atoms but infinite photons, reason why it is 
mandatory to set in numerical calculations a cutoff in. the 
photon Hilbert space. The convergence of our results is 
tested, checking their stability against small increases of 
this cutoff. 

Phase diagram.- As previously commented, two differ- 
ent phases, separated by A c , are found in Dicke model 
at zero temperature. The normal phase, where parity 
is a well-defined quantum number, and the superradiant 
phase characterized by a degenerated ground state and 
a spontaneous parity-symmetry breaking. A convenient 
method to see if this phenomenon is also present in ex- 
cited states is to analyze the difference 



Fig. [T] displays the dependence of these coefficients on 
N. Solid circles represent the numerical points corre- 
sponding to -Bjv; solid squares, the corresponding to An, 
and solid lines the fits to power laws N~ a . It is clearly 
seen that An — > —1 and Bn — > as N — > oo, and hence 
we conclude that E C (X) = E c (X,oo) = — J. It is worth 
mentioning that this value coincides with that recently 
obtained in the study of the connection of an excited- 
state quantum phase transition (ESQPT) with the de- 
velopment of quantum chaos and the critical decay of 
the survival probability 
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between the i-th excited states of both parity sectors 
II = ±1. If AEi is different from zero, the correspond- 
ing eigenstates have well-defined parity; if it is zero, they 
are degenerated and one can perform a rotation that mix 
both parity values. Results for N — 40 atoms are shown 
in Fig. [T] For A > A c there exists an abrupt change from 
AEi ~ to AEi > at a certain critical energy E C (X, N). 
A quantitative estimate of this energy can be obtained 
as the first eigenvalue B, for which AEi > fcern where 
fcerr is a given error bound. For all the results shown be- 
low we have set fc err = 10~ 6 ; similar ones are obtained 
with different bounds. Since the actual phase transition 
occurs in the thermodynamical limit, it is mandatory to 
infer how this critical line evolves as N — > oo. To do so, 
we assume that for each value of N the critical line obeys 
the linear law 



FIG. 1: (Color online). Intensity plot showing the decimal 
logarithm of the relative difference AE(E, X, N) in terms of 
E/J and A for a system with N = 40 atoms. Lighter regions 
correspond to the conserving-parity region, while darker re- 
gions represents the spontaneously broken-parity phase. The 
inset shows the finite-size scaling of parameters An and Bn 
inEq. ©. 

We can also use the mean-field approximation, that 
gives the exact ground-state energy in the thermody- 
namic limit, to estimate E C {X). Let us introduce for 
fi, v G M the coherent ansatz \fi, v) = ® \v), where 



Im) = (i + m 2 ) V j +|J,-J) 



\v)=e v /2 e va |0) , 



(4) 



correspond to the atomic and the photonic parts of the 
state, respectively. The resulting energy surface is 



E ° ) N =A N + B N X, A>A C . (3) E var Qi,u, \) = (u,p\H\u, 14 = 1*0 J 
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where the coefficients An and Bn are numerically de- 2 , \ ATT I ^ v 



termined by means of a least-squares fit. The inset of \ fi 2 + 1 
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It is plotted in Fig. [2j the upper panel shows the case 
with A = 0.25 (below A c ), and the lower panel, the case 
with A = 2.0 (above A c ); in both panels a number of 
level curves are drawn with solid lines . The geometry 
of this surface reveals that the level curve E = — J plays 
an especial role. For A < A c , it reduces to a single point 
at (/i, v) = (0,0), which is the absolute minimum of the 
energy surface; for A > A c , it changes abruptly to a non- 
analytic level curve containing a saddle point. Moreover, 
the shape of the energy surface is quite different depend- 
ing on whether E is below or above E — —J. In the 
former case, the energy surface exhibits two symmetric 
wells below E = — J, so that level curves are disjoint. 
On the contrary, for E > — J there is a single well with 
connected level curves for any value of A. This behavior 
supports that E C (X) = — J, as our previous numerical 
estimations for finite N. 




FIG. 2: (Color online). Contour plot of the energy surface 
E va , T (fi, v, A)/ J for two different values of A, one above and the 
other below the critical coupling A c . Upper panel, A = 0.25, 
and lower panel A = 2.0. Solid lines represent level curves. 

Dynamical symmetry breaking.- Baumann and co- 
workers [i| explored in real time the spontaneous parity 
breaking of the ground state at the superradiant phase 
transition, by measuring the behavior of ( J x ) as the cou- 
pling constant A increases in time and crosses the crit- 
ical point. Here, we follow an analogous procedure to 
study the different phases of the excitation spectrum 
when A > A c . We study the non-equilibrium dynam- 
ics and the relaxation to a steady state of certain physi- 
cal observables, like J x and q = (a* + a) /V%- They are 
physically measurable operators [j, [HI , which change the 
parity of the state on which they operate. Thus, they give 
rise to qualitatively different steady expectation values, 



depending on whether the energy of the non-equilibrium 
state is above or below E c . Although we only report 
results for J x , the behavior of q is completely similar. 
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FIG. 3: (Color online). Contour plot of the energy surface 
E(\i, A/)/J as function of Xi and A/. The critical energy line 
is placed at E/J = — 1. The darker region corresponds to en- 
ergies greater that the criitcal one, while the lighter zone cor- 
responds to lower energies. Solid line represent level curves. 

Let us take as our initial condition a symmetry- 
breaking ground state |^(0)) = |/Zj,z/j), where (/Uj,z/i) 
is the minimum of the energy surface corresponding to 
a coupling constant \ inside the superradiant phase. 
Then, we perform a diabatic change of A, i.e., a quench 
Aj — > A/, so that the ground state \(M,Ui) becomes a 
non-equilibrium state of H (A/). Its energy E(Xi, A/) = 
(fii, i>i\ H (A/) \fn, Vi) can be written as 

E(K,X f ) = ^ J (%) + 2J^—^ 

The contours of E(Xi,Xf) are shown in Fig. [3] It is 
clearly seen that choosing Xi and A / properly, one can ex- 
plore the different phases of the excited spectrum. In par- 
ticular, from any initial initial condition satisfying that 
Xi > 1.5A C , both phases can be reached by just quenching 
the system to different final coupling parameters A/. 

After performing the quench, we study the time evo- 
lution of (J x (t)) = (*(t)| J x |#(*)), where |$(t)) = 
e -iH{X f )t | vj> (0)) . If one expands the initial state in the 
eigenstate basis of H(Xf), denoted here as the 
expectation value of J x reads 

(■/*(*)) = E Cjde-W'- 111 * (Ej] J x IE,) , (7) 
being C t = (Ei\V(0)) . 
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FIG. 4: (Color online). Time evolution of the expectation 
value of J x , for two different quenches and N = 20. Green 
line corresponds to a time evolution with E < E c , and the 
red one to a evolution with E > E c . The inset displays the 
evolution at shorter times. 



Fig. @] displays, for a system with N = 20 atoms, 
the expected values (J x (t)) after applying two different 
quenches. In the first quench Xi = 1.41, A/ = 1.13, and 
the energy E(Xi,Xf)/J = —2.5 is well inside the parity- 
breaking phase. It is clear seen that (J x (t)) relaxes very 
quickly to a non-zero value. The same result is obtained 
in all the cases where the energy of the non-equilibrium 
state is E(Xi,Xf) < — J. Although in this particular 
case (Jx(t)) < 0, both positive and negative expectation 
values can be obtained. For every A^ > A c there are 
two degenerate ground states |±/Zj,±i/j), characterized 
by values of /i and v with opposite signs, that lead to 
different signs of (J x (t)). The starting point of the second 
quench is also A,; = 1.41, but the final coupling constant is 
reduced to A/ = 0.51. The energy of the non-equilibrium 
state, E(Xi, Xf) / J = 1.0, is now well above the critical 
energy. In this case, and in all that E > E c , we obtain 
(Jx(t)) = in the steady state. 

The physical explanation of this result is the follow- 
ing. For long-enough time evolutions, almost any ini- 
tial condition relaxes to a certain steady state, around 
which it fluctuates [3]. Moreover, when the energy 
eigenvalues are not degenerated, the expectation values 
of any observable O in the steady state are described 
by the diagonal approximation (0(t)) S ~ (O) D = 
J2i \Ci | 2 (Ei I O \Ei). On the contrary, if the energy eigen- 
values are degenerated, the diagonal approximation docs 
not hold, and thus it is possible that (0(t)) (0) D . 
These are precisely the cases that we have found in our 
model. In the preserving-parity phase [E > E c ) one can 
apply the diagonal approximation because the energy lev- 



els are not degenerated. As parity is a good quantum 
number in this case, (Ei\ J x \Ei) = for every energy 
level above E c . Therefore, whenever E > E c we find that 

(Mt)) ~* (J X ) D = Ei\Ci\ 2 (Ei\J*\Ei) = 0. On the 
contrary, in the broken-parity phase, the energy eigen- 
values are doubly degenerated with opposite parities, so 
that the diagonal approximation is not valid. Thus, one 
can find expectation values (J x (t)) in this phase. 

Consequently, the steady expected value of J x provides 
a neat signature of the two phases of the excited spectrum 
whenever A > A c . In fact, it acts like an order parameter 
of the ESPQT, as it is equal to zero if E > E c , and differ- 
ent from zero if E < E c . Therefore, it suffices to follow 
the long-time dynamics of a parity-changing operator to 
infer whether the energy of the initial state is above or 
below the critical energy. Furthermore, as this is already 
true for small values of N, finite precursors of this phase 
transition could be clearly observed in experiments. In 
particular, the setup used in Q is a good candidate for 
covering this aim. 

Conclusions.- We have studied the phase diagram of 
the Dicke model in terms of the coupling constant A 
and the energy E. Using numerical calculations and the 
mean-field approximation, we have found different phases 
in the excitation spectrum, separated by certain critical 
values E c and A c , where the latter also defines the critical 
point of the superradiant transition of the ground state. 
For A < A c we find a single phase where parity is a well 
defined quantum number. On the contrary, for A > A c 
there exists a critical energy E c — — J, such that below 
E c all the energy levels of the system are doubly degener- 
ated and composed of states with opposite parities. As a 
consequence, fluctuations can entail a spontaneous par- 
ity breaking — the system can evolve into a state without 
a definite value of the parity. In some sense the excited 
energy levels up to energy E c inherit the properties of 
the ground state in the superradiant phase. This fact 
leads to measurable dynamical consequences. Starting 
from a symmetry-breaking ground state in the superradi- 
ant phase and performing a quench to a non-equilibrium 
state, the relaxed expected value of certain observables, 
like J x or q, is different from zero only if the energy of 
the non-equilibrium state is below E c . This constitutes a 
new feature of the Dicke model, which could be observed 
in experiments similar to that of Ref. [H . We think that 
the results contained in this Letter might shed some light 
about the significance of the critical behavior of the Dicke 
model in real physical systems. 
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